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Abstract. - The Casimir interaction is omnipresent source of forces at small separations between 
bodies, which is difficult to change by varying external conditions. Here we show that graphene 
interacting with a metal can have the best known force contrast to the temperature and the Fermi 
level variations. In the distance range 50 — 300 nm the force is measurable and can vary a few 
times for graphene with a bandgap much larger than the temperature. In this distance range the 
main part of the force is due to the thermal fluctuations. We discuss also graphene on a dielectric 
membrane as a technologically robust configuration. 



Introduction. — The Casimir force [T] manifests it- 
self at short distances (< 1 (im) as a result of the elec- 
tromagnetic interaction between neutral bodies without 
permanent polarizations. For two ideally reflecting paral- 
lel plates separated by a distance a, this force is given by: 
Fq = (7T 2 / '240) (he /a 4 ). The universal character of the 
force stimulated active development of the field [5] with 
applications in physics, biology, and technology. 

The Lifshitz theory [3] gives the most detailed descrip- 
tion of the force. According to this theory, current fluctua- 
tions (quantum and classical) in the bodies are responsible 
for the force. Fluctuations in a wide range of frequencies 
give significant contribution to the force. For this rea- 
' son it is difficult to change the force at will as one has 
to modify the dielectric response of interacting materials 
in a wide range of frequencies. Hydrogen-switchable mir- 
rors did not show observable contrast to the Casimir force 
[4] . It was demonstrated that the force between indium tin 
oxide (ITO) and a gold surface is 50% smaller than it is be- 
tween two Au surfaces [5] . For the same material the best 
result was found for the phase-changing material (Ag-In- 
Sb-Te) with 20% difference between amorphous and crys- 
talline phases [6 . In situ modulation of the force between 
a gold sphere and a silicon membrane [7] was shown to 1% 
level when the carrier density was changed optically by 4 
orders of magnitude. 

The force measured in modern experiments is mainly 



the result of quantum fluctuations whilst the force due to 
classical fluctuations (thermal Casimir or Lifshitz force) 
was measured only recently between an ultracold atomic 
cloud and a sapphire substrate [Hj, and between two Au 
surfaces [9]. The thermal fluctuations dominate the force 
at large distances a > hc/T (Jcb — 1) where the force 
itself is extremely weak and approaches the Lifshitz limit 
[3]. Between two metals this limit is given by: 

T C(3) he 
F^^i a»A T = -, (1) 

where is the thermal wavelength and £(x) is the zeta- 
function. 

In this paper we show that significant variation (up to 
5 times) of the total Casimir force is possible for graphene 
with a bandgap 2A 3> T. The force changes in response to 
the variation of the Fermi level mainly due to the change 
of its thermal part. It can be realized at the distance range 
a = 50 — 300 nm, where the force is well measurable. 

Graphene, a single layer material with carbon atoms 
arranged in a honeycomb lattice, attracted enormous at- 
tention [lOjfTT] . Unusual electronic properties of graphene 
are due to massless relativistic dispersion of electrons at 
low energies [HJH2]. The Casimir /van der Waals interac- 
tion of graphene was mainly discussed at zero temperature 
[I3HT5] with the conclusion that the force due to graphene 
is weak in comparison with the interaction of bulk bodies. 
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An important development was made by Gomez-Santos 
[TB] at finite temperature. It was argued that at T = 
graphene is a critical system, with no characteristic length 
scale. At nonzero T this scale is given by the thermal 
length = hvp/T, where vp w 10 6 m/s is the Fermi 
velocity in graphene. It was found that in the long distance 
limit the force between two graphene sheets is given by the 
same Eq. (JTJ but this equation is true for much shorter 
distances a J> £t- At room temperature the scales fr 
and At are 25 nm and 7.6 (im respectively. This property 
makes the thermal Casimir force operative for separations 
in the 50 — 300 nm range, which are readily accessible 
using an atomic force microscope (AFM) or other force 
measuring techniques (see [2] for a review). 

Graphene is a promising material for the development 
of high-performance electronic devices |17j but pristine 
graphene is a semimetal with zero bandgap The 
major challenge of graphene electronics is to open an en- 
ergy bandgap [18]. As we will see later, the bandgap is 
also important for tailoring the Casimir force by electronic 
means. Significant progress has been made in this direc- 
tion. For instance, epitaxially grown graphene on SiC has 
a gap of 2 A « 0.26 eV [T5]. Opening a bandgap was also 
demonstrated by water adsorption [20] and patterned hy- 
drogen adsorption |21j . Graphene nanomesh proved to 
generate bandgaps with values depending on the mesh 
density 22-24 . Very recently, an efficient way to fabri- 
cate graphen nanomesh was developed |25) . In the present 
work we will assume the presence of a gap without speci- 
fying its origin. 

Graphene on a substrate. — We consider here the 
interaction between two plates 1 and 2 having dielectric 
functions £\(uj) and £2(w), respectively. In contrast with 
[16] graphene is not free standing but covers the plate 1. 
As we will see it has significant influence on the system. 
The case of suspended graphene is reproduced by taking 
£i(to) = 1. The Lifshitz formula [3] expresses the force 
between two parallel plates via their reflection coefficients. 
If graphene sheet has the two dimensional (2D) dynamical 
conductivity a, then the reflection coefficient of the plate 
with graphene (for p polarization) is given by [2"51l2"7] : 



k Q Ei — fci + (A-kcf/uS) k a ki 
fco£i + fci + (47TCr/w) fcofci 



(2) 



Here the normal components of the wave vectors in vac- 
uum and in the substrate are fco = \J u 2 j c 2 — q 2 and 
fci = yj £\ijj 2 1 1 c 2 — q 2 , respectively, where q is the wave 
vector along the plate. In the T = limit the graphene 
conductivity is a ~ e 2 /h for frequencies up to near UV 
[26,28,29. It means that the reflection coefficient gets 
only a small correction ~ a — e 2 /he = 1/137 due to the 
presence of graphene on the dielectric substrate. This ex- 
plains a weak force between two graphene sheets p~3lfl4] 
(2.6% of the force between ideal metals, ~ na). 

In this paper we neglect the effects due to a on the 
force. In this approximation the force between a sus- 



pended graphene sheet and any another material tends 
to zero at T = (negligible in comparison with the force 
between bulk materials). If graphene covers a substrate 
then the force difference AF = F g — Fj, is equally negligi- 
ble, where F g and Ft, are the the force with and without 
the graphene layer on the substrate, respectively. One 
can systematically neglect the effects ~ a in AF by tak- 
ing the non- retarded limit c — > oo. The possibility to use 
this limit was already indicated for two graphene sheets 
[IB] , Detailed calculation of the force between suspended 
graphene and Au [30] gave an independent proof of this 
approximation. Taking the limit c — > oo in the Lifshitz 
formula one finds the graphene contribution: 



/ oo 
n=0 £ 



R 



R 



— Rr> 



(3) 



where the integration variable in the physical terms is 
x = 2aq. Here R = t\T2 is the product of the reflec- 
tion coefficients for the body 1 (covered with graphene) 
and the body 2, and Rg = ror2, where r is the reflection 
coefficient of the body 1 without graphene. The reflection 
coefficients also have to be calculated in the non-retarded 
limit. The sum is taken over the imaginary Matsubara fre- 
quencies uj n = 2iirTn/h, which enter the dielectric func- 
tions in the reflection coefficients. Only p polarization 
contributes to AF since the s polarization vanishes in the 
non-retarded limit. It has to be stressed that c — > oo limit 
can only be applied to the force difference but not to F g 
or F\> separately. We keep the lower integration limit in 
J3J finite £„ = 2irTn(hc/2a)- 1 . Doing so we stay within 
acceptable uncertainty ~ a in AF. This definition is more 
convenient because convergence of AF is defined only by 
graphene but not high frequency transparency of the bulk 
bodies. 

To proceed further we need to know the dielectric func- 
tion of graphene. It is related to the dynamical conduc- 
tivity of the vacuum-graphene-dielectric system by the re- 
lation [3Tj : 



e{q,u) 



4na(q,uj) ( fc fc 



eifco + fci 



(4) 



Combining Eq. Q with Eq. @ one finds a simple ex- 
pression for the reflection coefficient of the body covered 
with graphene: 

e(q,u) 



Dielectric function of graphene. The dielectric 
function of graphene can be calculated using the random 
phase approximation (RPA). The RPA was used exten- 
sively for graphene in different situations (see the reviews 
[11][T2]). Specific to our case, we need to know this func- 
tion for imaginary frequencies at nonzero temperature for 
doped graphene with a nonzero gap. In the literature one 
can find s(q, of) only in different limiting cases. 
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For nonzero gap the electron energy in the valence 
(s = —1) or in the conduction (s = +1) band is Eg^ = 
syj (fiw^k) 2 + A 2 . The probability to find an electron 
(hole) with the energy E s k is given by the Fermi distri- 
bution / sk = [1 + e (s 8 fe-Ejr)/T]-i ) wnere e f i s tne Fermi 

level. In the RPA, the dielectric function of graphene can 
be expressed as e = l + v c (q)IL(q, u). Here v c = 2-ne 2 jnq is 
the 2D Coulomb interaction, k is defined by the environ- 
ment of the graphene layer (in our case 2k = £i(0) + 1), 
and H(q, oS) is the 2D polarizability given by the bare bub- 
ble diagram: 



n(g,w) 



-4E 



d 2 k 

W) 



rVi 



kk' 



/sk /s'k' 



(0) 



where k' = k + q, s,s' = ±1, and the vertex factor is 
given by 2Vg = 1 + (h 2 v 2 F k • k' + A 2 )/E sk E s , k , . The 
factor 4 at the front comes from two spins and two valleys 
degeneracy. 

In what follows we use the dimensionless variables: 



Q 



hvpq 
2T ' 



Z = 



2T' 



A T = 



T ' 



e F = 



Ep 
T 



(7) 



where C is the imaginary frequency. It is convenient to 
calculate the polarizability in the elliptic coordinates /i 
and v defined by the relations 

k = — (cosh /i — cos v) , k' = — (cosh /i + cos v) . (8) 

The notations £ = cosh /z and rj = cos v will also be used. 
Separating interband (k and k' in different bands) and 
intraband (k and k' in one band) transitions in ([B]) we can 
present the dielectric function as: 



e(g,<0 = l + — (A + /2) 



(9) 



where I\ and li are the contributions coming from inter- 
band and intraband transitions, respectively, and a g is the 
interaction constant in graphene. For / 12 one finds 



see eq. 



In eq. (T(| e 1;2 = \/Q 2 (£ TV) 2 + A T and the upper 
(lower) sign is related to index 1 (2). 

Typical values of q for the Casimir problem are ~ l/2a. 
Therefore, for the distances a > ^ of interest in this pa- 
per, the values of Q are always small, i. e. Q 1. In this 
limit eq. (|10|) can be simplified further. The parameter Qrj 
is always small but Q£ is not. In fact, the important values 
of £ in the integrals are large: £ ~ max(l/Q, At / Q) ■ Mak- 
ing the corresponding expansions and performing explicit 
integrations over v we find for a in the limit Q <C 1: 



7i = 7r<5 / <ie 



+ A T 



sinh e 



e 2 (Z 2 + e 2 ) cosh e + cosh cp ' 



(12) 
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v/Z 2 + g 2 - (A T Q/e) 2 
cosh e cosh 



(coshe + cosh 



(13) 



where we introduced a new integration variable e 
V Q 2 £ 2 + A T . Note that the intraband contribution d( 
inates the dielectric function in the Q <C 1 limit. 



The force. — In the large distance limit a 3> £t the 
dielectric function of graphene is significant (e - 1 » a) 
at frequencies fi-C S T, which are low for T ~ 300° K. For 
these frequencies most of dielectric materials have static 
permittivities and metals can be considered as perfect con- 
ductors. In such cases we can simplify the calculation of 
AF in eq. ([3]) taking the static permittivities £1,2(0) for 
bulk bodies (£2(0) — > 00 for metals) and keeping q and 
£ dependence only for the graphene dielectric function 
e(q,i(). It has to be mentioned that e(q,i() is essentially 
nonlocal. This nonlocality, however, is two dimensional, 
which simplifies the calculation of the Casimir force in 
comparison with the 3D case [32]. This is because there 
is only an in-plane wave vector. 

Consider first a gapless graphene. For A = the dielec- 
tric function at large distances a ^> £t follows from ([9]) 
and ([TBI 



£(<MC) = 1 + 



2a g G(e F , 0) 

Q 



Vz 2 



(14) 



The function G(ei?,0) here increases monotonously start- 
ing from 2 In 2 at = (fig. Q]i). In general G(x,y) is 
given by the expression 



G(x,y) 



dtt- 



1 + cosh t cosh x 
(coshi + coshx) 2 



(15) 



Let us stress that the characteristic frequency in the di- 
electric function (fl4| is £ ~ v F q as is expected from 
general consideration [TB]. For the Matsubara frequency 
cjo = *0 the dielectric function has a metallic character, i. 
e. e(q,i0) S> 1 and the reflection coefficient of the body 
covered with graphene approaches 1, i. e. r% — > 1. Al- 
ready, for n = 1 we have Z\ 3> Q and e(q,i(i) is strongly 
suppressed. For n ^ the reflection coefficient approaches 
the substrate value, n — ► ?*o- Let us stress that just one 
monolayer covering the substrate makes it perfectly re- 
flecting at low frequencies. 

For distances a > we can apply (fT4|) to calculate the 
force ([3]). The n — term dominates in AF. If the second 
body is a metal we can take R = 1 and i?o = ^0 > where tq 
has to be taken in the static limit. The force in this case 



is: 



AF(a,T) 



TC(3) 
87ra 3 



K(r ), a»^, (16) 
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h,2 = dfj, dv 



it 



Q(e 2 ±e 1 )(e 



4Z 2 +(e 2 ±ei) 



sinh 62 



± 



sinh 



cosh 62 + cosh ep cosh ei + cosh e f 



(10) 



where the function K(ro) describes the effect of the sub- 
strate on the force. This function is shown in figOJi and is 
expressed analytically as 



K(r ) 



1 



2C(3) 



dxx 



1 



>■() 



To 



(17) 



For suspended graphene ro = 0, and Eq. (|T6|) coincides 
with the Lifshitz force F^. Note that a metallic substrate 
for graphene will result in the zero force because K(ro) — > 
when ro — > 1. 

The effect of graphene will be appreciable if AF is mea- 
surable but also if AF is not negligible in comparison with 
the background force Fb. The force (fT6j) is maximal for 
free standing graphene when Fb = 0. This configuration 
is realizable in practice |34j and has significant interest. 
However, it can not always be practical due to the defor- 
mation induced by the force. A more stable configuration 
is graphene on a dielectric membrane of thickness h. For 
a membrane, the reflection coefficient is 



1 



-2qh 



ror, 



ro 



1 



-2qh • 



(18) 



where ro corresponds to the bulk material. For a thin 
membrane, h <c a, ro m becomes small and the background 
force Fb is much weaker than that for the thick substrate. 
For graphene-on-membrane the force in the long distance 
limit is also given by eq. (|T6|) but now the factor K de- 
pends slightly on the distance due to g-dependence of r 0m . 
The graphene-on-membrane configuration maximizes not 
only the absolute value of the force AF but also the rel- 
ative value AFjFb- This is an important practical obser- 
vation. 

Figure [Tb shows how the force approaches its limit value 
(fT6|) for free standing graphene, for graphene on 20 nm 
thick SiC>2 membrane, and for graphene on a thick Si02 
substrate. Numerical calculations were performed using 
the dielectric function ((9]) with 7i, 2 from (fl~0"|) without 
additional approximations. The continuous lines are for 
Ep = and the dashed lines are for Ep = 10 T. One can 
see that the force is not very sensitive to Ep. 

This is especially obvious in fig. [TJ; where the relative 
force (AF in respect to the background force Fb) is shown. 
This figure demonstrates significant dependence on tem- 
perature and shows that the relative force is considerably 
smaller for a thick substrate than for a thin membrane. 

Significant dependence on the Fermi level is desirable to 
change the force by electronic means. This can be realized 
if graphene has a non-zero gap. The material will change 
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Fig. 1: (Color online) (a) (top-right axes) Function G(x,0) 
that enters Eq. (|14J1 . (bottom- left axes) The factor K(ro) in 
Eq. (|16p as a function of the static reflection coefficient of the 
substrate supporting graphene. (b) The force ratio AF/Fl as 
a function of distance for free standing graphene, graphene-on- 
membrane, and graphene-on-substrate at T = 300 K. The solid 
lines are for cf = and the dashed lines are for ef — 10. (c) 
The relative force for graphene-on-membrane as a function of 
a for two different temperatures. The lowest dashed and solide 
curves are for graphene-on-substrate at T = 300 K. The curves 
for membrane in (b) and (c) were calculated for h = 20 nm 
and ro = 0.6. 

from insulating to conducting state in response to the po- 
sition of Ep. It has to influence the dielectric function and 
thus the force. The dielectric function of graphene with 
the gap 2A was calculated in [33 , on the real frequency 
axis at T — 0. Here we are using our result ©, (fTU|) for 
the dielectric function on the imaginary frequency axis at 
non-zero T. 

As in the case of gapless graphene the main contribution 
to the force at large distances comes from the n = term, 
which depends on the static dielectric function: 



In 

e(q,i0) = l+-^-G(e F ,A T ), 



(19) 



where the function G(ef, Ap) is given by eq. (TT51) . The 
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Fig. 2: (Color online) (a) The force as a function of the Fermi 
level for free standing graphene with At = 10. The thick 
line is for a = 50 nm and the thin solid line is for a = 100 
nm. The dashed, dash dotted, and dotted lines are the first 
three components for the a = 100 nm case, (b) The force for 
graphene-on-membrane (At = 10). The lines marked as 1, 2, 
and 3 correspond to a — 100, 200, and 300 nm, respectively, 
(c) The background force as a function of the distance in the 
units of the bare Casimir force Fc = iv 2 hc/240a 4 . 



gap gives significant effect for At 3> 1. If the Fermi level 
is in the middle of the gap, i. e. ep = 0, the function G is 
exponentially suppressed, i. e. G(0, At) ~ 2Atc~ At . In 
this case the effect of graphene on the force is small. When 
€f becomes comparable with At the dielectric function 
e(q, iO) is large and the effect of graphene is significant. In 
the long distance limit the force behavior is similar to eq. 
(Ull). 

Figure [5^ shows the force for suspended graphene with 
the gap At = 10 as a function of the Fermi level for 
a = 50 nm and 100 nm (solid curves). About ten terms 
are important in the sum ([3]); the first three terms for 
a = 100 nm are shown. Indeed, the n — term gives 
the main contribution. The finite value of the force at 
Ep — decreases as a and A increase. It is mainly due to 
intcrband transitions, which are not included in (|19j) . As 
expected, the force is small for Ep = and is on the level 
of Fl for the Fermi level Ep > A. Typically the force 
changes 3 — 5 times on the interval < Ep < A proving 
significant sensitivity to the Fermi level position. 

The force for graphene-on-membrane is shown in fig. ^jp. 
The behavior is similar to that for suspended graphene. 
However, in this case the force has to be compared with the 
background force for membrane shown in fig. [5J:. The lat- 
ter one was calculated using frequency dependent dielec- 
tric functions of SiC>2 and Au. The relative force AF/Ff, 



varies in the range 10 — 100%; it is small for short sep- 
aration and increases with a. The background force F^ 
can be reduced further by decreasing thickness and/or the 
permittivity of the membrane. 

Conclusions and discussion. In this paper we 
analyzed the Casimir interaction of a graphene-covered 
dielectric with a metal plate. The dielectric function of 
graphene was found at finite temperature for imaginary 
frequencies for the material with a finite bandgap and non- 
zero Fermi level. A simple expression ([3]) describes the 
graphene contribution to the force. We can conclude that 
for graphene with the gap 2A 3> T there is a strong depen- 
dence of the Casimir force on both the temperature and 
the Fermi level. This is realized at distances a Hvp/T 
when the main contribution to AF originates from ther- 
mal fluctuations. The predicted force is measurable with 
modern AFM instruments and can have significant tech- 
nological applications. Graphene-on-membrane interact- 
ing with a metal has special interest for practical appli- 
cations. This configuration combines mechanical strength 
with unique electronic properties of graphene. It allows 
tailoring of the Casimir force by electronic means. Ma- 
nipulations with the thermal force opens up completely 
new possibilities which, so far, seemed to have pure aca- 
demic interest for condensed matter. For example, it 
becomes possible to observe the nonequilibrium Casimir 
force [351 13S] between solid bodies at distances ~ 100 
nm. This possibility put the Casimir effect on the same 
ground as the short distance radiative heat transfer |37j . 
For all bulk materials the equilibrium component of the 
force at a ~ 100 nm is orders of magnitude larger than 
the nonequilibrium one. However, for suspended graphene 
or graphene on membrane interacting with a metal these 
components of the total force can be comparable. 
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